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For experiments with ions confined in a Paul trap, minimization of micromotion is often essential. In order 
to diagnose and compensate micromotion we have implemented a method that allows for finding the position 
of the radio-frequency (rf) null reliably and efficiently, in principle, without any variation of direct current 
(dc) voltages. We apply a trap modulation technique and focus-scanning imaging to extract three-dimensional 
ion positions for various rf drive powers and analyze the power dependence of the equilibrium position of the 
trapped ion. In contrast to commonly used methods, the search algorithm directly makes use of a physical effect 
as opposed to efficient numerical minimization in a high-dimensional parameter space. Using this method we 
achieve a compensation of the residual electric field that causes excess micromotion in the radial plane of a linear 
Paul trap down to 0.09 V/m. Additionally, the precise position determination of a single harmonically trapped 
ion employed here can also be utilized for the detection of small forces. This is demonstrated by determining 
light pressure forces with a precision of 135 yN. As the method is based on imaging only, it can be applied 
to several ions simultaneously and is independent of laser direction and thus well-suited to be used with, for 
example, surface-electrode traps. 

PACS numbers: 03.67.Ac, 37.10.Ty, 37.10.Vz 


I. INTRODUCTION 

Trapped laser-cooled ions are a prolific starting point for 
many experiments related to quantum information science 
iiiiia and precision spectroscopy, yielding some of the most 
accurate clock standards to date sa. Ions can be trapped 
for long times and laser-cooled, even to the motional ground 
state EHHl. They are one of the most promising candidates for 
quantum computation and for quantum simulations where 
the quantum mechanical properties of a system difficult to in¬ 
vestigate directly is simulated using the well-understood and 
-controlled quantum system of trapped ions sni. The size 
of this quantum system can be scaled up by entangling few- 
ion systems in complex segmented trapping architectures HB- 

[ni. 

The conceptual starting point for all this research is one 
or several ions at rest, more accurately, close to the motional 
ground state of an effective, approximately harmonic trapping 
potential ifTSl [T9l . The Laplace equation forbids the exis¬ 
tence of electrostatic potential minima in vacuum (Earnshaw 
theorem). This limitation is bypassed in radio-frequency (rf) 
quadrupole (Paul) traps, where the ions are confined by an 
inhomogeneous oscillating field ifTSl [T^ . and the oscillation 
energy of the forced oscillation becomes position dependent. 
Upon the presence of a cooling mechanism such as laser cool¬ 
ing, an ion comes to rest at the minimum of the oscillating 
field amplitude. Additional direct current (dc) fields from bi¬ 
ased electrodes or surface charges, e.g. , created by the load¬ 
ing process, can add forces which push the ion away from the 
rf null; in addition, a phase mismatch between rf electrodes 
might even prevent the existence of a time-independent rf null. 
In both cases, the ion will carry out a forced oscillation at the 
drive frequency, the so called micromotion, with an amplitude 
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depending on dc and phase mismatch. Micromotion plays an 
important role in excess heating of trapped ions, e.g., through 
unwanted Doppler shifts Il20l42^ and in conjunction with in¬ 
evitable anharmonicities or stochastic changes of the trapping 
potential ll24tl^ . 

Renewed interest in micromotion minimization comes from 
recent research with combined traps of neutral and charged 
atoms EMSl where excess micromotion is a detrimental 
source of neutral atom loss and has prevented collisional cool¬ 
ing of ions by a cold neutral atom cloud to the fundamental 
limits lf30l . Here, as an initial stage of combined traps, mi¬ 
cromotion is carefully compensated and the neutral atom loss 
can be used as a figure of merit for optimization tMl . Mi¬ 
cromotion can be overcome altogether by using a dipole trap 
for the ion. As the light interacts with the dipole moment 
of the ion instead of the charge, dipole traps are much softer 
and shallower and themselves require carefully balanced dc 
offset fields as a prerequisite ||32l. Another motivation for 
low micromotion comes from precision spectroscopy and fre¬ 
quency standards: micromotion sidebands spoil the accuracy 
of the determination of atomic resonance frequencies due to 
Doppler shifts, and thus micromotion minimization is an es¬ 
sential prerequisite for accurate ion trap frequency standards. 
One disadvantage, though not a matter of principle, is that ion 
clocks suffer in short term stability, as just one atom is probed. 
Massively parallel interrogation of many ions in common or 
separate traps, all micromotion compensated, is an important 
step to improve the short term stability of ion based frequency 
standards ll3^ . 

Micromotion minimization can thus be regarded as a com¬ 
mon crucial initial step and prerequisite for good trap perfor¬ 
mance. Micromotion causes diverse signatures, for which sev¬ 
eral minimization methods have been developed: 

1. The ion’s absorption and emission spectra are altered 
due to the periodic Doppler shifts caused by micromo¬ 
tion directly or by mechanisms conditioned on the pres- 
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ence of excess micromotion. 

Micromotion minimization is performed by exploiting 
the Doppler shift induced temporal correlation between 
the scattering rate and the rf trapping voltage ll3?l - [?5l 
and by measurement of micromotion sidebands in the 
absorption spectrum Ill5l - l39l or the emission spectrum 
Boll spaced by the trap drive frequency. 

For ion crystals the normal mode spectrum can be al¬ 
tered by excess micromotion and compensation of stray 
electric fields can be performed by minimization of 
these frequency shifts BD. 

2. The ion’s equilibrium position depends on the strength 
of the effective trapping potential if the dc fields are not 
nulled at the rf node. 

Excess micromotion is detected by monitoring the ion’s 
average position while changing the strength of the ef¬ 
fective trapping potential using either dc fields Il42l or 
changing the rf voltage amplitude llT5l . 

3. Parametric resonances of the ion’s motion can be ex¬ 
cited by a modulation of the rf trapping voltage, if the 
dc fields are not nulled at the rf node. 

The excitation of such a resonance causes a change in 
the ion’s scattering rate, which is then minimized ll43l - 

ia. 


sideband spectroscopy. A single position determination in 
two dimensions requires only milliseconds of measurement 
time, accumulating to a few hundred milliseconds of data ac¬ 
quisition for the complete minimization process. If position 
determination and the minimization are carried out in three 
dimensions, we need to use a tomographic method and the 
approach is slowed down by about an order of magnitude. On 
the other hand, even then, no constraints on the propagation 
direction of the laser limit the optimization procedure. As the 
position determination is based on imaging, the optimization 
might even be carried out for several ions in parallel, which is 
of particular interest for frequency standards. 

The method discussed here is fast, has a sensitivity compa¬ 
rable to that of other approaches and is yet general, simple, 
and economic. 


The paper is structured as follows: in Sec. we model 
the displacement of the ion’s equilibrium position. The ex¬ 
perimental setup is detailed in Sec. [TO The measurement pro¬ 
cedure and results are presented in Sec. IV As the position 
determination can be precise to the nanometer level, we use 
this for the measurement of small forces in the example of the 


light pressure in Sec. VI 


II. THEORY 


4. Micromotion causes an collisional transfer of kinetic 
ernergy from the ion’s micromotion to neutral atoms 
when trapped simultaneously. 

Micromotion minimization is performed by measure¬ 
ments of the loss rate of the atom trap EH. 

In this paper, we revisit the trap modulation technique, Q. 
For different rf levels, we determine the three-dimensional 
(3D) ion position from a tomographic imaging procedure. 
This position can be determined with an uncertainty far below 
the wavelength of light scattered off the ion for observation. 
From the ion positions measured at different rf powers, we 
can extrapolate the ion trajectory to infinite rf power, which 
ends at the rf null. With the knowledge gathered from a single 
trajectory, the ion can be moved to the rf null by changing the 
dc electric fields, which requires accurate a priori character¬ 
ization of the trapping fields. As an alternative, trajectories 
recorded at different dc settings readily yield this characteri¬ 
zation, together with the optimized compensation voltages. 

The time required for the minimization process can be cru¬ 
cial if the dc fields are fluctuating. This can happen as a conse¬ 
quence of stray charges and, usually to a lesser extent, by drift¬ 
ing fields from unstable voltage supplies ll44l |46l . For surface 
traps trapping times can be in the second or minute range ow¬ 
ing to smaller trap depths and might require frequent reload¬ 
ing and patch potentials may change on a short timescale. In 
such cases it is desirable to carry out the minimization fast and 
repeatedly to track or compensate time-dependent effects. 

The method discussed here relies on position determina¬ 
tion, which can be carried out with a high precision and does 
not require any narrow transition and ultra stable lasers for 


This section starts with an overview of the mathematical 
treatment of a single ion confined in a Paul trap in the pres¬ 
ence of an electric stray field, both by solving the respective 
Mathieu equations llT5l and by using the pseudopotential ap¬ 
proach lITSl . Readers familiar with this treatment might be 
referred to Eq. Q and Eq. (|^ as the main results and proceed 
from Eq. 0 for a discussion of the signature ion trajectories 
evolving from this models. 

An ion (mass m, charge Q) trapped in a Paul trap (operated 
at an rf voltage V cos fit and a dc voltage U) in the presence 
of an electric stray field obeys the inhomogeneous Mathieu 
equations ll35l . 

h + [ak+ 2gkCOs(flt)]^rk = ( 1 ) 

4 m 

where r = 'Y^rkk (k = x,y,z) is the ion position, k = 
X, y, z are the Cartesian unit vectors, ak and are the com¬ 
monly called trapping parameters, and the electric stray field 
Sstray = X] Astray,is approximated to be static and uniform 
over the small volume occupied by the ion trajectory. 

This stray field might be a real physical field (e.g., originat¬ 
ing from charged insulators near the trap), but also effects lift¬ 
ing ideal symmetry and displacing the dc origin with respect 
to the rf (e.g., misaligned trap electrodes or machining imper¬ 
fections) as well as other homogeneous static forces (e.g., the 
light pressure force exerted by a cooling laser) can effectively 
be treated as contributions to the stray field. 

Solving Eq. Q using the adiabatic approximation 
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(\ak\,q^ ^ 1) yields the approximate motion of the ion lf35l 

rk{t) « [re,k + fo.fc cos(wfcf + ipk)] (l + y cos{Vtt)^ (2) 

= [ro,fe cos{uJkt + ipk)] (l + y cos(f^f)^ 

+rE,k + cos(nf), (3) 

where Tq and ipk depend on the initial conditions and ujk = 
^y/ak + ql/2. Eq. § describes a motion characterized by 
two frequency components: the secular motion component, 
oscillating at frequency uJk, and the micromotion component, 
oscillating at the much higher frequency and much lower 
amplitude qkro,k/2. The secular motion can be reduced by 
cooling mechanisms as, e.g. , laser cooling. This reduces the 
micromotion amplitude as well, as it is proportional to the 
secular motion. The last two terms in Eq. 0 describe the 
effects of the additional electric field. The field Sstray shifts 
the average position of the ion out of the rf potential node to 
the position 


re 



^stray,fc 



k, 


(4) 


at which the rf electric field causes oscillations with ampli¬ 
tudes qkre,k/2 - the excess micromotion. This motion is a 
driven motion and thus cannot be significantly reduced by 
cooling techniques. 

In the pseudo- or ponderomotive potential approach im, 
the motion of the ion due to the rf electric field is averaged 
over one period of micromotion and the so-called rf pseudopo¬ 
tential represents the kinetic energy due to micromotion. 


^:f{r) 


m 

Y 


k 


2 

^rf./c ^ 2 ■ 


(5) 


The Wrf^/c are called the rf trap frequencies and represent the rf 
contribution to the secular frequencies ujk- 

The ion can then be treated as being confined in an effective 
potential that is the sum of the rf pseudopotential and all 
dc potential contributions [including dc electrodes (</)dc,i) as 
well as stray potentials ((/^stray)], 


= (/>rf(r) +Q'^ ^dc,^(r) + Q^stmy(r), (6) 

i 

where —V^stiay(t’) = Astray Here the are taken to be 
ideal potentials in the aforementioned sense: contributions 
causing a shift of the dc origin relative to the rf node (that 
is, contributions linear in the ion’s position) are absorbed into 

Astray ■ 

Near the center of the confining potential (e.g. , along the 
axis of a linear trap) the sum of dc potentials can be approxi¬ 
mated as a quadrupole potential: 

Qj2(IJdcAr) = y (7) 

i k 

Similarly to the Wrf fc, the uJdc,k are the dc trap frequencies at 
a specific setting of all dc trapping voltages (wf^, ^ = H^afe/4 




FIG. 1. (Color online) Left: In the presence of a stray field Estray, 
for example, an electric stray field from misaligned electrodes or the 
effective field of a light pressure force, the ion is pushed out of the 
minimum of the effective rf potential (thin blue line) to ro where it 
experiences a nonvanishing rf field and thus is subject to micromo¬ 
tion. Right: When the applied rf power is changed, the curvature of 
the effective rf potential and the equilibrium position of the ion 
change. As the stray fields in the left and right plots are of identical 
strength, the equilibrium condition as given in Eq. is fulfilled for 
identical slopes (dashed gray lines) of the effective rf potential at the 
ion’s position. 


in the case of an ideal linear trap). Eor the sake of simplicity 
the principal axes of the rf and dc potentials have been chosen 
here to be parallel and equal to the Cartesian axes. 

The equilibrium position of the ion is then given by zero net 
force 


0 = -V^eff(r) 
= -V 

Q 


y ^2 + ‘^dc.fc) rl + QAtrayir) 

k 

^stray,fc 


( 8 ) 


^0 = — V 2 , 2 

k + ‘^dc.fc 


Astray,/c 


k 


_ Q _ 

m t ^ 




(9) 


which is equivalent to the average displacement from 
Eq. 0, as ujI = Here P^f/Po has been in¬ 

troduced as a convenient scaling factor. Prf is the applied rf 
power, and uJrfg^k the rf trap frequencies at an applied rf power 
of Pq (P (X V'^ (x q\)- The shift of the equilibrium position 
out of the rf potential minimum depends on the stray field and 
the applied rf power (see Eig.[2. 

As an intuitive point of view, it is also useful to combine the 
effects of all stray fields (real as well as effective fields) as a 
translation ^(estray) of the dc potential’s origin explicitly and 
rewrite the effective potential as 

<t>aii{r) = y ^ (10) 

+ y ^dc.k [rk — <^fc(£stray)] • (H) 

k 


In a ID trap, the electric stray field would shift the dc po¬ 
tential by 


(5 


Q Astray 


mu!, 


dc 


( 12 ) 
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and the ion’s position in the trap by 


frequencies are given by 


ro 


Q 

m 


Po 


^ stray 


2 ■ 

dc 


(13) 


2 _ Q" 2 

2m2f22“'=- 


(15) 


As the applied rf power is increased the ion’s position 
will then approach the minimum of the rf pseudopotential 
(ro = 0). By monitoring ro as a function of the rf power 
an extrapolation of the rf node in the limit of Prf —oo is 
possible. 

Instructive cases of 2D potential configurations displaying 
characteristic features of the ion’s movement are shown in 
Fig. 0 The electric stray potential has been included as a 
shift of the dc potential’s origin [see Eq. (fTT)]. Rf potential 
isolines are drawn in yellow; dc potential isolines, in blue. 
Filled colored circles indicate the respective potential’s ori¬ 
gin. The thick black line marks the equilibrium positions of 
the ion while tuning the applied rf power from infinity to zero. 

Fig. displays a configuration of confining rf and dc po¬ 
tentials with degenerate trap frequencies (oj^f ^ j, > 0, 
X — ‘^dc y > 0)’ realized in the central hori¬ 

zontal plane of a ring trap with a positive dc bias on the ring 
electrode. Due to radial symmetry, the ion’s equilibrium po¬ 
sitions will move in a straight line between the rf and the dc 
node as the rf power is varied. 

'When radial symmetry is broken (e.g. , ^ y) 

the ion’s trace changes to a hyperbola (Fig.|^). The dashed 
line indicates the full hyperbola and its asymptotes. As the 
rf potential still features rotational symmetry, the hyperbola’s 
asymptotes are parallel to the dc potential’s principal axes. 

For now we have assumed the rf and dc potentials to be 
confining in all directions. In a typical linear Paul trap, the 
dc potential will be used to create dc confinement in the axial 
direction. Since the dc potential has to obey Faplace’s equa¬ 
tion [A(/)dc('t') = 0], it must be repulsive in at least one of 
the radial directions. Fig.|^ shows the same configuration as 
Fig- 03, but with an unstable dc trapping in the x direction 
(<^dc X < 0)’ resembling the potentials in the transverse plane 
of a linear Paul trap. The extreme values of the two poten¬ 
tials do not reside on the same branch of the hyperbola but 
are seperated by its pole. When decreasing the rf power and 
such lowering the rf confinement, the ion’s equilibrium posi¬ 
tion will move away from both, the rf minimum and the dc 
saddle point, finally escaping from the trap when the rf con¬ 
finement becomes weaker than the dc repulsion. 

In the most general case (Fig. |^) also the rf potential’s ra¬ 
dial symmetry is broken ^ y) and the directions of 

the dc potential’s principal axes are rotated with respect to the 
rf potential. This causes a shearing and rotation of the hyper¬ 
bolic trace, featuring nonperpendicular asymptotes. 

A particularly interesting consequence of nondegenerate ra¬ 
dial rf trap frequencies is the existence of a nonvanishing axial 
rf component. In the quadrupole approximation the rf potenial 
can be written as 

^'^ak{V)cos{nt)rl, (14) 

k 

where the depend on the actual geometry of the trap and 
the rf voltage V. Using Eq. s and qk = the rf trap 


Since the rf potential fullfills Faplace’s equation 
t) = 0 Vf, it follows that 

A(/)[“'(r, t) =2 cos(Uf) ^ ak = 0 

k 

^^afc = 0. (16) 

k 

A vanishing axial rf component requires ctz = 0, which 
immediately yields the degeneracy of the radial rf trap fre¬ 
quencies, 

CXx = ~OLy => Wrf,a: “ ^rf,y (^"^^ 

The above considerations allow us to detect excess micro¬ 
motion by monitoring the ion’s position while varying the ap¬ 
plied rf power. Using Eq. (j^ (adjusted for a rotation of the 
dc potential’s principal axes), in principle, the measurement 
of a single trace of the ion’s movement during an rf variation 
without any variation of dc voltages is sufficient to extrapolate 
the location of the rf node and therefore the location of zero 
excess micromotion. 

In practice, insufficient knowledge of the trapping parame¬ 
ters (e.g. , relative rotation of the rf and dc potential’s princi¬ 
pal axes, electric field generated by compensation electrodes) 
might require the combination of variations of the rf power 
and several dc voltage settings. An example of such a com¬ 
bination is shown in Fig. The rf and dc potential config¬ 
uration is the same as in Fig. |^, with three additional set¬ 
tings for compensation voltages, shifting the origin of the 
original dc potential (yellow lines) to (1.0, —0.5) (blue lines), 
(—0.5, —0.5) (green lines) and (—0.5,1.0) (red lines). The 
dc potential’s origins are marked by the respectively colored 
dots. The ion’s traces (thick lines) for all of these settings con¬ 
verge to the rf node in the limit of infinite rf power. A parallel 
analysis of the measurements for different dc configurations 
then yields the position of the rf null as well as the required dc 
compensation voltages. 


III. EXPERIMENTAL SETUP 

A schematic of the region of the microstructured segmented 
linear Paul trap used in the experiments presented here is 
shown in Fig. The trap electrodes have a separation of 
500 pm in the y direction and 350 pm in the x direction. The 
dc electrodes are divided into segments of 250 pm width (la¬ 
beled dcup,i and dcLow,; in Fig. |^. The full trap features 33 
pairs of dc electrodes composing a wide section (described 
above) and a narrow section (electrode width, 100 pm; y sep¬ 
aration, 250 pm) that are connected by a tapered transfer sec¬ 
tion. The trap and the supporting experimental setup have 
been described in detail in ll^ l47l l48l . 
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FIG. 2. (Color online) Two-dimensional potential configurations of the dc (blue isolines) and rf (pseudo-)potentials (thick yellow isolines). 
The potentials’ origins are indicated by respectively colored filled circles. The thick black line represents the ion’s trace as the strength of 
the rf potential is tuned from oo to 0. The dashed gray line shows the full hyperbola and its asymptotes describing the ion’s movement, 
(a) A configuration of degenerate trap frequencies y > 0, = ^dc,y > 0). Due to radial symmetry, the ion’s equilibrium 

positions will move in a straight line, (b) For broken dc radial symmetry (culcx < ‘^dc.y) the ion’s trace changes to a hyperbola, (c) A 
potential configuration with an unstable dc trapping in the x direction (cj^cx < 0) as realized in the transverse plane of a linear Paul trap, (d) 
Additionally broken rf radial symmetry > a;^,y) and a rotation of the dc potential’s principal axes causes a shearing and rotation of the 
hyperbolic trace, featuring nonperpendicular asymptotes. 



FIG. 3. (Color online) A combination of different dc potential 
shifts for a potential configuration as shown in Fig.|^. Rf (pseudo- 
)potentials isolines are shown in gray. The ion’s traces (thick lines) 
are shown for the original dc potential (yellow lines) and the shifts 
of its origin to (1.0, —0.5) (blue lines), (—0.5, —0.5) (green lines) 
and (—0.5,1.0) (red lines). Dashed lines show the respective full hy¬ 
perbolas; respectively colored filled circles the dc potentials’ origins. 
All traces converge to the rf node (filled black circle) in the limit of 
infinite rf power. Near the rf node, the set of all trajectories resem¬ 
bles a star, as expected by intuition. The directions of the hyperbola’s 
asymptotes are unchanged by the shifts. 


For this work we trapped single ions, which are 

produced by a two-photon ionization process from neutral yt¬ 
terbium vapor emitted from the atom oven 14^ . The ion is 
Doppler cooled at the S 1/2 ~ Pi /2 transition. The ionization 
laser (399 nm), the cooling laser (369 nm), and two additional 
repumpers (638 nm and 935 nm) are overlapped outside the 
trap BtII and propagate through the center of the trapping re¬ 
gion, enclosing angles of approximately 35° with the x direc¬ 
tion and 89° with the y direction (see Fig.|^). 

For basic trapping the central electrode pair (dcup/Low,o) is 
set to —IV, providing axial confinement, whereas all other 


dc electrodes are set to ground. The rf electrodes are driven 
through a helical resonator with a voltage of 210 Vpp at a 
frequency of 13.2 MHz. The secular frequencies of an ion 
trapped with this configuration are Wa,/ = 27r x 1.054 MHz, 
LOy! = 27r X 1.112 MHz and = 27r x 199 kHz. The princi¬ 
pal radial axes of the trapping potential {x' and y') are rotated 
by about 64° with respect to the x and y axes. 

Static and approximately uniform fields for micromotion 
minimization are created using three compensation voltages 
Uc (seeFig.gr). Uci is added as an antisymmetrical bias volt¬ 
age to all upper (h-Uci/2) and lower (-Uci/2) dc segments, 
causing a field roughly along the x + y direction near the 
trap axis. A second compensation voltage, Uc 2 , is applied 
to the aperture of the nearby atom oven. The aperture extends 
mainly in the y-z plane and its size largely exceeds the dimen¬ 
sions of the central trap slit, creating a field mainly along the 
X direction. A bend of the aperture along the y direction near 
the trapping site causes a non-negligible z component of the 
compensation field created by Uc 2 (see Fig.|^). Therefore a 
counteracting field is created by adding a calibrated bias volt¬ 
age Ub(Uc 2 ) antisymmetrically to all outer left (dcup/Low,i, 
z < — 2) and outer right (z > 2) dc electrodes (colored or¬ 
ange), which cancels the unwanted z component of the aper¬ 
ture’s field. The third compensation voltage Uc 3 is applied 
antisymmetrically to the same outer electrodes if control of 
the ion’s axial position is desired. 

The ion is detected by imaging fluorescence of the cooling 
transition along the —x direction on an EMCCD (Andor iXon 
X3 Blue) (see Fig. |^): The light gathering objective (numer¬ 
ical aperture of 0.4) ISOl is mounted on a motorized traveling 
stage and allows a positioning of the focus of the imaging sys¬ 
tem with a resolution of 8 nm along the x axis. The objective 
collects the ion’s fluorescence with a magnification of approx¬ 
imately 12.5 onto the chip of the EMCCD, consisting of 512 
X 512 pixels with an edge length of 16 pm x 16 pm. 
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FIG. 4. (Color online) Experimental setup (not to scale): (a) Top 
view (along the —y direction) including laser beams crossing the 
trap (enclosing angles of approximately 35° with the x direction) 
and showing the bend of the oven aperture. The oven aperture is not 
shown in (c). (b) Side view (along the z direction) showing the up¬ 
per and lower electrode layers (dc and rf) and the Yb oven aperture 
used as a compensation electrode (voltage Uc 2 ). The second com¬ 
pensation voltage is applied by antisymmetrically biasing all upper 
and lower dc electrodes (voltage Uci). Additionally a schematic 
of the focus-scanning 3D determination of the ion’s position [see 
Sec. l |IV A| l] is shown. The light gathering objective is mounted on 
a motorized traveling stage, which is used to sweep the focal plane 
(vertical dashed black lines) while imaging a series intensity distri¬ 
butions (solid blue lines). The ion’s x coordinate is determined from 
the fit of a Gaussian beam waist to the spot sizes at the different fo¬ 
cus settings (thick dashed red lines). The y and z coordinates are 
taken as the center of a 2D Gaussian fitted to the intensity distribu¬ 
tion imaged closest to the x coordinate.; (c) Frontal view (along the 
X direction) of the upper and lower dc electrodes (colored) and rf 
electrodes (semitransparent gray). 


IV. METHOD AND RESULTS 
A. Eocus-scanning position determination 

Since our experiment features imaging of the ion only from 
a single direction, we determine the full 3D position of the 
ion in our trap during the rf and dc variation by employing 
focus-scanning imaging (see Fig. 0’)- For every position de¬ 
termination, we take several images of the ion with the fo¬ 
cus swept over ±15 pm. For every single image a 2D Gaus¬ 
sian is fitted to the intensity distribution, yielding (y, z) center 
coordinates of the distribution as well as the Gaussian spot 
size. The x coordinate of the ion is then deduced by fitting 
the spot size of a Gaussian beam to the spot sizes obtained 
from the single images (see Fig. [^. Then the center coordi¬ 
nates of the intensity distribution imaged with the focus set¬ 
ting closest to the ion’s x coordinate are taken as the y and 
z coordinates of the ion. Using this procedure, one is able 
to determine the average ion position with a precision that is 
significantly below the wavelength of the fluorescence light. 
In principle, the precision scales as the inverse square root of 
the observed fluorescence and is limited only by its signal to 
noise ratio and the mechanical stability of the imaging sys¬ 
tem. For an average number of about 7000 photons per image 
(observed scattering rates of about 200 kHz to 400 kHz) and 
images taken at 25 focal positions, we routinely achieve a pre¬ 
cision of (TyjiJz) = (98 nm, 39 nm, 35 nm). 


B. Micromotion minimization in two dimensions 

To minimize micromotion in the radial plane at a spe¬ 
cific axial position we choose a set of voltage settings 
{{Uci,Uc 2 )i] and apply a set of rf powers {Pifj}- For all 
combinations {i,j) a focus-scanning position determination is 
carried out, yielding the ion’s equilibrium positions 
To suppress possible systematic effects of a quasicontinuous 
change in the rf power, we iterate over the different rf powers 
not in ascending or descending order, but over some permuta¬ 
tion of the chosen set. 

We then model our effective trapping potential as a super¬ 
position of the rf pseudopotential, the unshifted dc potential, 
and two additional linear potentials as created by our compen¬ 
sation electrodes, 

- rrf) 

+ '^ir - r^fRa,,M^2R~l{r - r^) 

±Scl(t^Cl,i - Uc\)y'^Racii'f' - f’rf) 

±SC2(C^C2,i — Uc2)x^Racii.'’'' ~ ^rf), (18) 

where R^ is a rotation by a, M ^2 and A<U 2 are the diago- 

^rf ^dc 

nal matrices of the squared rf and dc trap frequencies ^ 
and Wjj, ^), and sci and sc 2 are geometrical factors depending 
on the compensation electrode shape and placement. Uc\ and 
Uc 2 are the compensation voltages that cancel the stray po¬ 
tentials and shift the origin of the dc potential to the position 
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(a) 


No. of photons 



z (px) z (px) z (px) 


(b) 



FIG. 5. (Color online) Focus-scanning 3D determination of the ion’s 
position, (a) Fitting t2D Gaussians to the intensity distributions taken 
at different focus displacements yields the center coordinates (y,z) of 
the distribution as well as the spot size. Exemplary fits are given for 
focus displacements of ±13.5 pm and 0 pm. White lines indicate 
the distributions isolines at one and two FWHM. (b) The ion’s x 
coordinate is deduced by fitting the spot size of a Gaussian beam 
waist to the spot sizes obtained from the single images. Error bars 
for large focus displacements are overestimated as the assumption of 
a gaussian beam only holds near to the focus. The center coordinates 
of the intensity distribution imaged with the focus setting closest to 
the ion’s x coordinate are taken as the ion’s y and 2 coordinates. 


of the rf null (r^). The trap frequencies are determined from 
independent measurements by excitation of the ion’s motion 
with a voltage modulation of one dc segment for various set¬ 
tings of the applied rf power Pjf. All remaining parameters are 
determined by fitting the equilibrium positions obtained from 
Eq. ( [T^ to the 

Fig. shows the results of an optimization procedure with 
12 settings for the compensation voltages {Uci , Uqi) and vari¬ 
ation of the applied rf power over 12.5 dB in 25 steps (varying 
the rf voltage from 420 Vpp to 99 Vpp). The circles represent 
the ion’s positions tq i j with color-coded rf power, the dashed 
line represents the fit, and the position of the rf node Trf is 
marked by the central black circle. After applying the opti¬ 
mized compensation voltages as yielded by the fit, we perform 
a similar measurement but without any variation of the com¬ 
pensation voltages. While lowering the rf confinement, we de- 
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EIG. 6. (Color online) A 2D micromotion minimization performed in 
the radial plane perpendicular to the trap axis. The x axis is along the 
line of sight of our imaging optics. Optimization was carried out us¬ 
ing 12 settings of compensation voltages (Uci, Uci) and 25 rf power 
settings. Colored filled circles represent the ion’s positions with their 
respective rf powers, and example error bars indicating ±cra; and 
±2 ctj, are given for four points in the upper-right quadrant. Eit re¬ 
sults of the ion’s trajectories are indicated by dashed lines and the rf 
node is marked by the filled balck circle. The residual electric stray 
field uncertainty along the principal radial axes after the optimization 
is determined to (Asstray.a:', Acstray.y') = (0.09 V/m, 1.09 V/m) at 
respective trap frequencies of = 27r x 132.5 kHz and uiyi = 
2ti X 387.0 kHz, where the principal axes of the trapping potential 
{x', y') are rotated by 64.8° counterclockwise with respect to the x 
and y axes. 


tect no change in the ion’s average position before it escapes 
the trap (at an rf voltage of about 62 Vpp), where the escape 
happens more rapidly than the time resolution of our position 
determination. Using Eq. ([^ and the weighted standard devi¬ 
ation of the ion positions along the principal radial axes of our 
trapping potential, we derive a residual electric stray field un¬ 
certainty of (Asstray,a:', Asstray.y') = (0.09 V/m, 1.09 V/m) 
at respective trap frequencies of Wa,/ = 27r x 132.5 kHz and 
ujy> = 27r X 387.0 kHz, where the principal axes of the trap¬ 
ping potential {x', y') are rotated by 64.8° counterclockwise 
with respect to the x and y axes. 


Rf power dependence of the trap structure 

During the measurements the rf drive power is varied, so 
the energy deposition at the trap structure will be changed. 
The accompanying change in temperature might affect the 
trap geometry due to thermal expansion and might result in 
an rf power dependence of the rf node. The rf node yielded 
by the trajectory analysis might then be different from the rf 
node present during normal trap operation. 

If the time between changes in the rf power is short com- 
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pared to the time constant for temperature changes in the trap 
structure, the rapid changes in the rf drive will average to an 
effective rf power throughout the measurement. Thus the tem¬ 
perature of the system does not change significantly if the ef¬ 
fective power during the measurement is matched to the con¬ 
stant rf power during normal operation. This matching can 
always be achieved by the choice of rf powers that are sam¬ 
pled during the minimization process. 

In our particular setup the trap is composed of gold-plated 
aluminium oxide glued and wire-bonded to an actively cooled 
aluminium oxide holder BtI . Given the low thermal resistiv¬ 
ity of our trap-holder setup and its good thermal coupling to 
the cooling system, we estimate the time constant for temper¬ 
ature changes of the trap-holder system to be large compared 
to the time between changes in the rf power for our measure¬ 
ments (a few tens of milliseconds). Also, our comparatively 
low maximal rf trapping voltage (420 Vpp) (and, thus, low 
heating power) limits the system to only small temperature 
changes. In addition, the order in which measurements at dif¬ 
ferent rf powers are carried out adds to the averaging effect. 
Thus, the conditions for an effective constant rf drive power 
are fulfilled in the experiments presented here. 

A non-negligible change in the trap geometry caused by the 
changes in the rf power would manifest itself in a discrepancy 
between model and measurement data. A change in the trap 
geometry would yield both a change in the rf trap frequen¬ 
cies [Eq. Q] not compatible with the changes predicted 
by the model and a change in the dc trap frequencies Wdc.fc 
[Eq. (|7]l] that are assumed constant. These effects would alter 
the ion trajectory as a function of the rf power [Eq. (j^], which 
is not observed in the experiment (Fig. |^. 

In general, even if a trap structure exhibits a time constant 
for rf-drive-induced temperature changes that is comparable 
to the time required to take a single image of the ion, and 
if thermal deformation of the trap geometry is of concern, 
the latter can be avoided. The single-image exposure can 
always be shortened such that the averaging condition men¬ 
tioned above is fulhlled. The whole measurement sequence 
will then be executed repeatedly and the single-image posi¬ 
tion determination is performed by averaging over repetitions 
of the measurement sequence, with only a minor effect of in¬ 
creased readout noise. 

If constant rf power is mandatory, the method presented 
here can also be carried out by a variation of a dc voltage 
bias applied to the rf or dc electrodes instead of the rf power 
variation ll42ll . Under variation of a dc bias symmetrically ap¬ 
plied to the rf or dc electrodes, the rf trap frequencies will 
remain constant, while the dc trap frequencies change. As in 
the case of an rf power variation the ion’s equilibrium position 
will follow hyperbolic trajectories [Eq. (|^]. 

C. Micromotion minimization in three dimensions 

By construction, an ideal linear trap features no axial rf 
electric held component and therefore no axial micromotion. 
For all realizations with hnite-size electrodes, particularly for 
segmented traps, regions of non-negligible axial micromo¬ 


tion cannot be avoided, although recent results have demon¬ 
strated trap designs with small residual axial rf held compo¬ 
nents 1331 muni EH. 

The focus of segmented traps, such as the one used here, of¬ 
ten lies in the storage and manipulation of ion crystals. Such 
operations, e.g., shuttling, splitting, and merging, make use of 
an extended region along the almost-micromotion-free sym¬ 
metry axis of the linear trap. In these cases a space curve of 
positions of minimized radial micromotion along the axial di¬ 
rection may be used to optimize the operation of the trap. In 
a similar sense, when experimenting with extended linear ion 
crystals oriented along this axis, hnding a single point of van¬ 
ishing axial micromotion is not meaningful. Any axial offset 
from this ideal position, as long as it is much smaller than the 
extent of the crystal will have a small effect. So the level of 
precision aimed for in the determination of the axial position 
will be of the order of a micrometer or slightly less. 

If the axial positions of trapped ions are not constrained 
by other experimental requirements, minimization of axial as 
well as radial micromotion can be performed. A 3D micromo¬ 
tion minimization is demonstrated qualitatively by extending 
the procedure for 2D minimizations by a variation of the axial 
compensation voltage Uc 3 - The model of the effective trap¬ 
ping potential given in Eq. ( [T^ is complemented by adding 
the respective terms for axial components of the rf pseudopo¬ 
tential, dc potential, and linear potential created by Ua- Fig. [7] 
demonstrates a 3D optimization procedure performed in a re¬ 
gion of our trap suffering from signihcant axial rf electric field 
strength. The minimization was performed with 4 different 
settings for the radial compensation voltages {Uci, Uqi) and 
3 different settings for the axial compensation voltage Ucs, 
yielding 12 ion trajectories of varying rf powers. Filled cir¬ 
cles represent the average ion positions, with their colors in¬ 
dicating the respective rf powers. A fit of the trajectories from 
the 3D effective trapping potential are given by dashed lines. 
As described in Sec. (j^ all trajectories approach the rf node 
(filled black circle) for increasing rf field strength. 


V. PERFORMANCE 

In EH a comparison of several micromotion minimization 
approaches is given in terms of various figures of merit, such 
as the residual electric stray field uncertainty, the micromo¬ 
tion amplitude, and the average kinetic energy of micromo¬ 
tion, which relates to the second-order Doppler shift. For a 
specific experiment, one of these numbers might be more suit¬ 
able or relevant than the others. The micromotion amplitude 
and the average kinetic energy depend, to different extents, on 
the ion mass and on the drive and secular frequencies. These 
numbers vary substantially throughout the ion trapping com¬ 
munity, therefore the residual electric stray held uncertainty 
Ae, from which the other quantities mentioned above can be 
derived, is the relevant quantity that can easily be compared 
among experiments. In general, by increasing the rf drive 
power, one can, within the constraints given by the experi¬ 
ment, reduce the effect of micromotion. 

In table (|I| the outcomes of several minimization methods 
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FIG. 7. (Color online) Results of a 3D optimization procedure per¬ 
formed in a trap region with non-negligible axial rf field strength. 
Minimization was performed using 4 (Uci,Uc 2 ) x 3 (Ucs) settings 
of compensation voltages. Filled circles represent the average ion 
positions obtained from focus-scanning imaging, with their colors 
indicating the respective rf powers. Example error bars indicating 
itajj and ±2(jy_2 are given for four points. A fit of the trajectories is 
indicated by dashed lines. The filled black circle shows the position 
of the rf node. 


are listed. The final level of precision of micromotion min¬ 
imization can be scaled in each method by changing the av¬ 
eraging time for data taking, as long as the stray fields to be 
compensated are stable. Ultimately, as our minimization ap¬ 
proach relies on imaging, the position determination and, in 
turn, the determination of the rf null location are shot noise 
limited. Thus another question to ask is the sensitivity of a 
minimization method, that is, on what timescale what level of 
precision is obtained. Here, almost no information was found 
in the literature. 

The observed fluorescence rate is limited by the S' 1/2 — ^i /2 
transition of '^^Yb'^and its linewidth T = 27r x 19.2 MHz, 
the solid angle covered by our light gathering optics (4 % 
of 47 r) and its transmission, and the quantum efficiency of 
our camera (55 %). With these constraints and the time re¬ 
quired for the translation of our optics (about 400 s), record¬ 
ing a data set such as shown in Fig. takes about 620 s and 
directly yields all compensation voltages. We achieve a sen¬ 
sitivity of the residual electric field in the radial directions of 
= (2.3 V/mViU, 27.1 V/mViU). 


VI. POSITION DETERMINATION APPLICATION: LIGHT 
PRESSURE MEASUREMENT 


Another application of precise position measurements is the 
determination of small forces, where trapped ions have been 
proven to be excellent probes 1531 iMl- We demonstrate this 
with the example of the light pressure force Il55l acting on a 
laser-cooled ion in a Paul trap. 

When micromotion is minimized using any scheme that uti¬ 


lizes a signal from an incoherent scattering process, the light 
pressure force due to the fluorescence-inducing laser is com¬ 
pensated automatically. As a consequence, the settings ob¬ 
tained from the minimization are only valid at the scattering 
rate used during the process and will be different if the inter¬ 
action is changed (e.g. , change in laser intensity or use of a 
different transition). This is especially true for coherent ma¬ 
nipulation of the ion, as there is no light pressure force present. 

If minimized micromotion is required for situations without 
photon scattering, it is neccessary to perform the minimization 
for several scattering rates and extrapolate the minimization 
settings to the zero-scattering level. Given sufficient knowl¬ 
edge about the dc potential, one can monitor the ion’s position 
as a function of the scattering rates and counteract solely the 
light-pressure-induced shift. 

The force corresponding to the time-averaged momentum 
transfer of absorbed laser photons (or light pressure force) F\p 
is given by 


Fip = {pph)^=TflK, (19) 


with rate of absorption F, Planck constant H, and laser wave 
vector K. 

This force can effectively be treated as originating from a 
stray potential and adds as Fip k /w to the right-hand-side of 
Eq. Q. In analogy with Eq. (Q, the light pressure force shifts 
the ion’s average position by 

rip = -Mj/Fip , (20) 

m 


with being the inverse matrix of the squared secular fre¬ 
quencies. In the frame of the potential’s principal axes M ~2 


diagonal and Eq. (20i reduces to f ip = ^ X) [see 


Eq. (|4|]. 

When the light field interacting with the cooling transition 
is detuned by 5 from resonance, the scattering rate F and 
hence the rate of absorption of the ion are given by 


■So _ 7/2 

1 + so 1 + (2(5/702 


( 21 ) 


with the on-resonance saturation parameter sq = 2172 / 7 ^, the 
Rabi frequency H of the interaction, the natural line width of 
the transition 7 , and the saturation broadened line width 7 ' = 

7v/l + So- 

The expected light pressure force of a laser saturating 
(so = 1) the Si /2 — Pi /2 transition of red detuned 

by 7/2 is ji^ipl « 35 zN. For the parameters given as set¬ 
ting (a) in Fig. this force would cause a shift of rip « 
(4 nm, 1 nm, —5^0 nm), which is of the order of the level of 
precision of the position determination reported above for the 
z component and significantly below for the x and y compo¬ 
nents. To be able to observe the light pressure shift, averag¬ 
ing over several position determinations is a neccessity in our 
present setup. Furthermore, the mechanical stability of our 
imaging system requires differential measurements to cancel 
slow movements of the objective mount causing a virtual shift 
of the ion’s position. 












10 


TABLE I. Comparison of micromotion minimization methods in terms of the residual electric stray field uncertainty Ae. The first six rows of 
data are an extension of the values given in ED; the last row states the results obtained in this work. If different residual electric stray fields 
for different trapping directions are reported, only the lowest value per reference is given. 


Ref. 

Method 

Ae(V/m) 

HD 

Photon-correlation spectroscopy 

0.9 

ID/f38l/07l 

Micromotional sideband spectroscopy 

7 / 1 / 0.4 

da 

Ion-cavity emmission spectroscopy 

1.8 

(451/(41 

Parametric excitation of secular motion 

6/0.4 

ED 

Neutral atom loss 

0.02 

m 

Monitor displacement 

< 11.8 

This work 

Trajectory analysis 

0.09 


We measure the shift of the ion’s center coordinates in the 
z-y plane as a function of the change in the observed scatter¬ 
ing rate for different sets of trap frequencies. The ion’s scat¬ 
tering rate is varied by tuning the intensity of the 369 nm light 
field using an acusto-optic modulator in 15 steps. Exposure 
times are adjusted such that the amount of collected photons 
for each intensity setting is approximately constant to avoid 
systematic effects in the position determination. We obtain 
the shifts of the ion’s central positions {Ay, Az) with respect 
to the lowest intensity setting and average these shifts over 500 
repetitions of the intensity variation for each trapping config¬ 
uration, which improves the precision to the nanometer range. 
The change in observed scattering rate ATobs is determined 
by extracting the number of collected fluorescence photons 
from the ion images and referencing to the lowest intensity 
setting as well. Using a part of the image with negligible con¬ 
tributions from the ion’s intensity distribution we monitor and 
correct for the background. 

Fig. shows the shifts of the ion’s average position in y 
(blue symbols) and z direction (yellow symbols) versus the 
change in observed scattering rate for three sets of trap fre¬ 
quencies: (a) strong confinement in all three directions, (b) 
weakened axial confinement, and (c) weakened radial confine¬ 
ment. Detailed trap frequencies and corresponding symbols 
are given in the legend in Fig. The inset shows a zoom-in 
where data from the weak axial confinement have been left out 
for clarity. The linear dependence of the ion’s position shift on 
the observed scattering rate and thus the light pressure force 
given by Eq. ( |20| ) and Eq. ([T^ is clearly visible. 

We fit a linear model, Arip^fc = ak ■ ATobs + Pk {k = 
y, z), to the data, where Uk — 7 ^ 7 is a 

scaling factor accounting for the overall efficiency of our 
imaging system (T = 7 • Fobs)- When we weaken the 
axial confinement of the ion but keep a comparable radial 
confinement [going from case (a) to (b)], the movement in 
z direction increases from = —0.19(3) nm/kHz to 
^ = —4.66(4) nm/kHz whereas the shift in the y direction 
= —0.01(1) nm/kHz) changes only negligibly. Keep¬ 
ing the axial confinement but lowering the radial confinement 
[going from case (a) to (c)], we observe no change in the 
shift along the z direction but an increased shift along y of 
= 0.144(13) nm/kHz, which agrees with the approxi¬ 
mate orientation of the trapping potential’s principal axes. 



ATobs (kHz) 


FIG. 8. (Color online) Shifts of the ion’s y (blue symbols) and « (yel¬ 
low symbols) position as a function of the change in observed scat¬ 
tering rate for three trap configurations. Corresponding symbols and 
trap frequencies (in kHz) are specified in the legend. Set (a) gives 
strong confinement in all three directions; set (b) features weakened 
axial confinement resulting in an increased shift along the z direc¬ 
tion. Set (c) is weakened in terms of radial confinement and causes a 
substantial shift along the y direction. Inset: Zoom-in on the region 
of small shifts. Straight lines are fits of the linear dependence of the 
position shift on the change of the light pressure force (see text). 


Since M ~2 is not diagonal in our frame of observation, 
the ion’s displacements cannot in general be associated with 
only one of the trap frequencies, and as the principal axes of 
the dc and rf effective trapping potentials differ, changing the 
radial confinement also changes the orientation of the radial 
trapping axes. For the shifts of the ion position along the z 
direction, however, a- = 7 —%, as the axis of our linear 
trap is aligned in parallel with z. So the ratio of for two 
different axial confinements should be equal to the squared 
inverse ratio of the trap frequencies Ratios obtained from 
the measurements displayed in Fig. and the respective trap 
frequency ratios for the axial position shifts are listed in table 
yielding a good agreement. 

The best sensitivity and uncertainty are achieved for 
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TABLE II. Comparison of the ratios of as obtained from our fits 
and the squared inverse ratios of the trap frequencies lJz (oiz oc 
for three trapping configurations (a), (b) and (c) (see text and Fig.[^. 



(a),(b) 

(c),(b) 

(c),(a) 


0.046(3) 

0.046(3) 

0.9947(15) 

ja^z'* 

0.041(5) 

0.044(5) 

0.92(16) 


lightweight particles and low trap frequencies along the di¬ 
rection of the scattering force [see Eq. As the sensitiv¬ 

ity depends on position accuracy, the measurement can again 
be optimized for either precision or speed. For our setup we 
report a sensitivity for the measurement of the light pressure 
force along the z direction of SF,p = 633 yN/ v/Hz (total expo¬ 
sure of a differential measurement: 38 ms), with our minimal 
absolute uncertainty being — 135 yN. 

VII. CONCLUSION AND OUTLOOK 

In our present setup, fluorescence is detected normal to the 
plane of the trap chip and, thus, perpendicular to the soft trap¬ 
ping axis along which confinement is almost purely realized 
by dc fields, labeled as the axial direction of our linear seg¬ 
mented trap (compare Fig.|^. In an ideal linear trap, the rf ef¬ 
fective potential and thus micromotion depend exclusively on 
the radial coordinates, perpendicular to the axial direction. In 
this case position changes upon rf variation will occur only in 
radial directions and thus partially in the direction of our line 
of sight. In real traps, there might be a residual axial rf elec¬ 
tric field, which is minimized by careful design and typically 
much smaller than the radial rf field strength. In segmented 
traps, the intention is usually to exploit the entire axial span 
of the linear trap by shuttling, etc., and thus axial micromotion 
might be present but potentially cannot be minimized. Micro¬ 
motion minimization always requires a determination of the 
radial ion displacement or motion, in our case also along our 
line of sight. We realize this with focus-scanning imaging as 
detailed above. In our present setup this is one of the main 
limitations, in both precision and speed. As our detection of 
excess micromotion is based on a shift of the ion’s average po¬ 
sition, the method as presented is insensitive to micromotion 
induced by a phase difference between the rf trap electrodes 
or an rf pick up of the dc electrodes that is phase-shifted by 
nearby filter circuits. These effects can be avoided by design 
to a great extent, especially in surface traps, and are often of 
negligible size. 

Detection along the axial direction of the trap would al¬ 
low us to obtain the relevant displacements without the time- 
consuming physical translation of massive light gathering op¬ 


tics and reduces the number of required images considerably: 
right now, we take 25 focal steps and need to include addi¬ 
tional times to wait for the completion of a translation action. 
Additionally, omitting a translation stage further improves the 
mechanical stability and thus the precision of the position de¬ 
termination. Taking into account translation and settling times 
as well as the reduction in the number of images, the speed 
of acquisition would increase by a factor of approximately 70, 
meaning that the entire minimization as shown in Fig. would 
be finished in two dimensions in roughly 9 s and the sensitiv¬ 
ity of this method would increase to about 0.3 V/m-\/Hz for 
a 2D minimization. Speed not only is relevant to keep the 
time for calibration, etc., low and leave most for the exper¬ 
iment of interest, but also allows to compensate more often 
and characterize and counteract time-dependent drifts of the 
power supplies and patch potentials. In principle, when the 
effects of dc field changes of each electrode are well charac¬ 
terized, a single trajectory will suffice, and this might further 
improve the performance of the method in terms of speed. 

Utilizing the same technique, precise position determina¬ 
tion of a harmonically trapped ion may also be used for the 
simple detection of small forces. We have demonstrated the 
detection of forces on the yoctonewton (10“^^ N) scale in the 
example of the light pressure force exerted on the ion by the 
cooling laser. Depending on the experiment, precise knowl¬ 
edge of the light pressure force might be desirable to counter 
its effects and maintain minimal micromotion also under ex¬ 
perimental conditions without incoherent scattering. 

In conclusion, we have presented a method, which allows 
us to compensate micromotion due to dc stray fields in three 
dimensions using focus-scanning imaging as well as the de¬ 
tection of rather small forces. The method is general and does 
not require spectroscopy of a narrow transition or a specific 
orientation of laser beams, which is of particular interest for 
surface-electrode traps. If imaging is carried out along the 
symmetry axis of a linear trap and the optimization is carried 
out in the transverse direction only, no additional components 
and no focus-scanning imaging are required and the optimiza¬ 
tion can be carried out within seconds. 
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